Measurement of Photon Statistics with Live Photoreceptor Cells 
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We analyzed the electrophysiological response of an isolated rod photoreceptor of Xenopus laevis under 
stimulation by coherent and pseudo-thermal light sources. Using the suction electrode technique for 
single cell recordings and a fiber optics setup for light delivery allowed measurements of the major 
statistical characteristics of the rod response. The results indicate differences in average responses of 
rod cells to coherent and pseudo-thermal light of the same intensity and also differences in signal-to- 
noise ratios and second order intensity correlation functions. These findings should be relevant for 
interdisciplinary studies seeking applications of quantum optics in biology. 
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Eyes of living organisms represent advanced light har- 
vesting systems, developed through hundreds of millions 
years of evolution. Some of their features are compa- 
rable or even superior to existing man-made photodc- 
tcction devices. For example, rod photoreceptor cells 
of the retina, which are responsible for night vision and 
form the focus of the present study, represent miniatur- 
ized photodetectors containing a photosensitive element 
(rhodopsin pigment) along with a "built-in" chemical 
power supply (ATP produced by mitochondria) [l| . They 
have sensitivity down to single-photon level, and demon- 
strate a remarkable low-noise operation 0]. Under- 
standing such properties of nature-given photodetectors 
stimulates considerable interest in interfacing them with 
sources of non-classical light, such as light with a "fixed" 
number of photons, and "squeezed" light Q. Implemen- 
tation of such "bio-quantum" interfaces would contribute 
to a better understanding of the performance of the visual 
system near the threshold [3] , allow a reference- free cali- 
bration of its quantum efficiency [f| , along with a direct 
detection of entangled multi- photon states 0. Study of 
the statistics of rod responses to excitations by classical 
light sources allows careful characterization of rods, and 
open possibilities for future studies of non-classical light. 

Vertebrate retinal photoreceptor rod cells, further re- 
ferred to as rods, are densely packed in the retina laying 
the inner surface of an eyeball Individual rods con- 
sist of two distinctive functional regions: extended rod- 
like outer segment (OS) which is filled with rhodopsin 
photopigment molecules, and shorter rounded inner seg- 
ment (IS) which contains various components of cell ma- 
chinery. In the dark, the electrochemical gradient across 
the rod cell membrane causes a continuous flux of ions 
(Na + , K + ,Ca 2+ ) in and out of the cell through aque- 
ous pores in the membrane, referred to as ion channels. 
A photon, impinging on the OS, causes isomcrization of 
a rhodopsin molecule, which initiates a chain of intra- 
cellular reactions, referred to as phototransduction cas- 
cade. A fraction of ion channels closes, thus resulting 
in hyperpolarisation of the cell. Illumination with a rel- 
atively bright flash closes all light sensitive channels in 



the OS, and bleaches the rod. A recovery process regu- 
lates re-opening of the ion channels to the resting state. 
Accurate measurement of the ion current flux through 
the rod membrane is possible with the suction-electrode 
technique 0. In this case, the OS of the rod is drawn in 
a tight fitting glass microelectrode filled with the phys- 
iological solution. The rod functions in a way similar 
to in-vivo system, however the ion current through the 
membrane is now re-directed to a low noise amplifier. 

In earlier experiments isolated rods were mainly stim- 
ulated by multi- mode thermal sources (lamps, LEDs), 
whose statistics, however, was not within the scope of in- 



terest of those works [8Hll|. Influence of well controlled 



light statistics on the response of the entire visual system 
was studied only in behavioral experiments with human 
subjects 0, EH. However, noise properties of individual 
rods are difficult to access in this case, as several hun- 
dred of rods are simultaneously illuminated, and their 
collective response undergoes several intermediate stages 
of the visual processing. In the present study, influence 
of controllable photon fluctuations on the response of an 
isolated rod is analyzed for the first time. Using well 
characterized coherent and pseudo-thermal light sources, 
together with a fiber optics setup for light delivery al- 
lows us to determine major statistical characteristics of 
the rod in the whole dynamic range of its response. 

The probability distribution of photon numbers is 
given by [l3j : 
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where m is the average number of photons. Based on 
earlier findings, it is further assumed that each photon 
isomerizes only a single rhodopsin molecule fl4| . In this 
case the probability of isomcrization of n molecules by 
m impinging photons is given by the binomial distribu- 
tion: P(n\m) = ( , ™)t7™(1 — rj) m ~ n , where 77 is defined 
as an overall quantum efficiency of the photodetection 
process [l3| . Then, the probability of isomerization of n 
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molecules by the fluctuating light source is given by: 
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It is further assumed, that in the absence of saturation 
subsequent electrophysiological responses to individual 
isomerizations are additive, and have a standard Gaus- 
sian shape Q. In this case, the probability of observing 
the amplitude of the photocurrent A in response to n 
isomerizations is given by: 
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where Aq is the average, and a a is the standard deviation 
of the photocurrent amplitude in response to a single 
isomerization, and ctd is the standard deviation of the 
photocurrent in darkness. Hence the final probability 
of observing the photocurrent amplitude A is given by 
averaging over Pi (n) : 



P(A) = ^2,P{A\n)Pi(n) 



(4) 



In order to account for the saturation of the photocur- 
rent, P{A) is truncated at the saturation amplitude 
As, so that the probability of observing the amplitude 
A > As just adds up and results in observing As'- 



Ps(A) 



'P(A), if A < A s 

^-lAlo P ( A ) dA ^ XA = A s 
0, if A > A s 



(5) 



Eq. (5) allows calculation of any statistical mo- 
ment of the photocurrent amplitude, given by A k = 
J A=Q A k Ps(A)dA. The relation between the average 

A and the variance VarA = A 2 — A 2 , can be conve- 
niently characterized by a dimcnsionless signal-to-noise 
ratio SNR = A/^Var(A). 

Investigation of joint statistics of the rod response with 
some reference detector, measured in a Hanbury-Brown 
and Twiss (HBT) type experiment 15, HH, allows finding 
the second order intensity correlation function of the light 
field g( 2 \ The value of g^ may evidence non-classical 
properties of light [T^]. From the joint pulse-to-pulse data 
of photocurrent amplitudes of the rod and photocounts 
of a conventional avalanche photodiode (APD), g^> is 
found as: 
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where K = rjDrn is the average number of APD photo- 
counts, i]d is the APD quantum efficiency, and P(K) is 
the probability distribution of APD photocounts, which 
is similar to (1). 
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FIG. 1: (color online) (a) Optical layout. A laser beam is 
chopped by a shutter, attenuated by a neutral density filter 
(NDF), and split by a 50/50 beamsplitter (BS). Both of the 
beams are coupled into single mode optical fibers (SMF) with 
lenses (L). One of the fibers is connected to an avalanche 
photodiode (APD) whilst the other stimulates the rod. Inset 
shows the pseudo-thermal light source: L is a lens, GD is a 
rotating ground glass disc, D is a diaphragm, (b) Microscope 
image (top view) of the fiber taper and a suction micropipette 
with the constrained rod cell. 



In the experiment (see Fig. la) the beam of a 50 mW 
continuous wave frequency-doubled Nd:YAG laser at 532 
nm (Photop, DPL-2050) was chopped by an optical shut- 
ter (Melles Griot) in 30 ms duration pulses with a repe- 
tition rate of 10-15 seconds. The pseudo-thermal source, 
shown in the inset of Fig. la, was realized by focusing 
the beam of the same laser with a lens (/=100 mm) into 
a rotating ground glass disc (GD, 1500 grit), followed 
by a 0.5 mm iris diaphragm (D), placed at a distance 
1200 mm [l7j. The diameter of the diaphragm was set 
to be much smaller than the size of a single scattered 
speckle. Then, either the laser beam (for the experiment 
with the coherent source) or a speckle (for the experi- 
ment with the pseudo-thermal source) was attenuated by 
a variable neutral density filter (NDF) , and directed to a 
HBT setup, consisting of a non-polarizing 50/50 beam- 
splitter (BS) and two single mode optical fibers (SM 460- 
HP) with equalized coupling efficiencies preceded by as- 
pheric lenses (/=6.24 mm). One of the fibers was used 
to stimulate the cell, whilst the other one was connected 
to an actively quenched APD (Perkin Elmer, SPCM- 
AQRH-15-FC). The number of photons used to stimu- 
late the cell was calculated from the number of APD 
photocounts in each pulse, assuming the APD quantum 
efficiency t}d =50% jl8|. The statistics of coherent and 
pseudo-thermal sources are checked in independent mea- 
surements by the APD (see supplementary materials). 

Methods of cells preparation, electrophysiology record- 
ings, and light coupling were similar to the ones we de- 
scribed previously [19j . In brief, dark adapted adult male 
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frogs Xenopus laevis were used for the experiment. All 
procedures with the animals were carried out in accor- 
dance with IACUC regulations. Shortly after the sacri- 
fice, the retinas were taken from the eyes under IR-light 
(850 nm) and sliced with a blade to release single rods. 
The rods were loaded into a Ringer solution filled cham- 
ber of an inverted microscope (Leica, DMI 3000), which 
was placed in a light tight Faraday cage. The micro- 
scope was equipped with an IR LED (920 nm), and a 
CCD camera (Leica, DFX 380) to allow observation of 
the preparation without bleaching the rods. A glass mi- 
cropipette, pulled to the diameter of ~5/um and bended 
by about 45° was used for conventional suction electrode 
recordings Q- OS of intact rods were captured in the 
pipette by gentle air suction, see Fig. lb. The photocur- 
rent pulses were amplified and digitized by a patch-clamp 
amplifier (Heka, EPC 10), the data were stored by the 
computer. The experiment was synchronized by an ex- 
ternal pulse generator, which triggered data acquisition 
by the rod and the APD, and conducted at room tem- 
perature. Each cell was used for continuous recordings 
during 90-130 min; 1-2 cells per preparation were used. 
The consistency of the response during the experiment 
was checked prior to acquisition of every point by excit- 
ing the saturation amplitude As with a dim torch Q. 

A tapered single mode optical fiber (Nanonics, working 
distance 25-30 /im, spot size 4 fim) , mounted on a motor- 
ized stage (Sutter, resolution 0.3 /xm), was used for light 
delivery to the rod (see Fig. lb). The fiber tip and the 
rod were carefully positioned at a distance of 25-30 nm, 
so that the stimulus were sent along the rod long axis. 
Such an "axial configuration" provides a nearly optimal 
match between the light spot and the rod profile, long 
absorption path in the rod, and an invariance of the rod 
response to light polarization [lj| . 

Actual waveforms of the rod photocurrent in response 
to stimulation by different numbers of impinging pho- 
tons, produced by the coherent source, are shown in the 
inset of Fig. 2. In addition, a histogram of rod responses 
to dim coherent light flashes is shown in Fig. 6 of supple- 
mentary materials. The dependence clearly indicates the 
quantization of the response, which supports the notion 
that rods function as photon number resolving detectors 
[H, in accordance with the suggested model. The de- 
pendencies of the normalized average photocurrent am- 
plitudes, on the average numbers of photons for coherent 
and pseudo-thermal sources are shown in Fig. 2. Ampli- 
tudes were normalized by the saturation value As, and 
the numbers of photons were normalized by mo, the value 
which initiated a response of a half saturation amplitude, 
which, for the majority of studied cells, was in the range 
of 550-2,500 photons/pulse. Lines represent numerical 
calculations of the averages, with the probability distri- 
bution given by (5). The dependencies for the major- 
ity of investigated cells yielded parameters in the range: 
A s e[18,26] pA, A) e[0.5,0.85] pA, a A e[0. 25,0.6] pA, 
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FIG. 2: (color online) Dependence of the normalized average 
amplitude of the rod photocurrent on the average number of 
impinging photons for coherent (red squares, solid line) and 
pseudo-thermal (black circles, dashed line) sources. Each data 
point is an average response to 80-100 light pulses. 6 different 
cells, taken from 6 different animals were used. Lines show 
numerical calculations, and standard deviation is shown by 
error bar. Inset shows raw waveforms of the rod photocurrent 
in response to coherent pulses of different average intensity. 

a D £=[0.1-0.2] pA, and r) e[0.007,0.01]. Note that r) ac- 
counts for an overall efficiency of the photodctection in a 
given experiment. 

The dependencies in Fig. 2 demonstrate that satura- 
tion of the average amplitude for the coherent source is 
significantly steeper than for the pseudo-thermal one. A 
clear explanation for the difference can be found from the 
experimental histograms of the normalized amplitudes, 
which are shown for both of the studied light sources in 
Fig. 3. At relatively low number of impinging photons 
rod operates in a linear regime, and the statistics of its 
responses displays the statistical features of the incoming 
light (Fig. 3 a,c). However, when the rod is stimulated by 
bright flashes, its saturation causes damping of response 
fluctuations and modifies the light statistics (Fig. 3 b,d). 
For the coherent source photon numbers are well local- 
ized around their average value (see (1)), and the ma- 
jority of response amplitudes to bright flashes are close 
to the saturation amplitude As (Fig. 3b). In contrast, for 
the pseudo-thermal source with the same average photon 
number, the probability of emission of pulses with rela- 
tively low photon numbers remains non-negligible (see 
(1)). Hence the amplitude probability distribution gets 
significantly broader than in the case of the coherent 
source (Fig. 3b), resulting in smoother saturation of the 
average response, as shown above (see Fig.2). 

Note that the capability of the entire visual system of 
living human subjects to distinguish between the light 
sources of different statistics was shown in Our ex- 
periments suggest that responses of single rods are at 
the basis of this phenomena. Our results also imply that 
the light source characteristics should be carefully con- 
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FIG. 3: (color online) Probability histograms of normalized 
amplitudes measured for coherent (a,b) and pseudo-thermal 
(c,d) sources at different values of the average, shown with 
the corresponding SNR values. Each histogram was acquired 
for 80-100 light pulses. Standard deviation is shown by error 
bar. 



sidered for faithful comparison of experiments on visual 
transduction in rods. 

The effect of saturation is evident from the experimen- 
tal dependence of the SNR on the number of impinging 
photons, shown along with numerical calculations using 
(5), in Fig. 4a. The raw data for the measurements of 
SNR, obtained for 6 different cells, are shown in sup- 
plementary materials (see Fig. 7). In the linear regime 
of the rod response, SNR grows linearly for the coherent 
source, and remains constant for the pseudo-thermal one. 
At the same time, it demonstrates a steeper growth in the 
saturation region for both of the sources. Note that sat- 
uration of rod should be considered in experiments with 
non-classical light. In particular, saturation might lead 
to experimental artifacts, such as observation of sub-shot 
noise fluctuations (" squeezing" ) in experiments with clas- 
sical light sources 

The experimental dependencies of the correlation func- 

(2) 

tion <7q on the number of photons for coherent and 
pseudo-thermal light sources, along with numerical cal- 
culations, are shown in Fig. 4b. The result for the coher- 

(2) 

ent source yields g y a '= 1.01 ±0.01 and does not demon- 
strate any dependence on the number of impinging pho- 

(2) 

tons. In contrast, for the pseudo-thermal source, g$ 
decreases with the increasing number of photons. Thus 
" excess" correlations of the thermal source are gradually 
neglected due to the saturation of the rod response. As 
an additional check, in the experiment with the pseudo- 
thermal source and for the linear regime of the rod re- 

(2) 

sponse, g)- was calculated between the photocurrent am- 
plitude and the number of APD counts for consecutive 

(2) 

light pulses. One observed a drop from g^ =1.94±0.07, 
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FIG. 4: (color online) Dependence of the SNR (a) and <?q 
(b) on the average number of impinging photons for coherent 
(red squares, solid line) and pseudo-thermal source (black cir- 
cles, dashed line). Results are averaged similarly as in Fig.2. 
Standard deviation is shown by error bar, which can not be 
seen for the coherent source in (b). Lines show numerical 
calculations. 



pulses are uncorrelated. Thus it is shown that rods al- 

(2) 

low adequate measurement of g^ similar to conventional 
man-made detectors, and could be readily used for char- 
acterization of quantum light. 

In conclusion, being inspired by the ultimate charac- 
teristics of rod photoreceptors, we carefully investigated 
the impact of photon fluctuations of various classical 
light sources on their response. The experimental re- 
sults revealed capabilities of isolated rods in measure- 
ment of photon statistics. It is of future interest to in- 
vestigate rods interfaces with nonclassical light, in partic- 
ular with correlated two-photon light and intense multi- 
photon twin-beam states. The developed approach can 
be also used in interfacing rods with realistic sources of 
intensity fluctuations, for example blinking star lights, 
when observed through turbulent atmosphere (20| . 
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SUPPLEMENTARY MATERIALS 
Source calibration 

Coherent and pseudo-thermal sources were character- 
ized by measuring the dependence of the variance of the 
APD photocounts on their average. The probability dis- 
tribution of APD photocounts is equivalent to (1), with 
m — >• K , m —> K. For the coherent source Vax(K) = K, 

and for the thermal source Var(AT) = K + K* 0. The 
experimental results, shown in Fig. 5, were found to be 
in good agreement with the theoretical model, validating 
the experimental approach. 
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FIG. 5: (color online) Dependence of the variance of the 
APD photocounts on their average for the coherent (full red 
squares), and the pseudo-thermal (open black circles) light 
sources. Lines show theoretical calculations. Gray shadow 
region shows the range of values which corresponded to the 
excitation of a half-saturating amplitude of the cell photocur- 
rent, denoted as mo in the text and in Figs. 2, and 4. 



Photon number resolution by the rod 



Conventionally single photon responses of isolated 
rods manifest themselves by quantized amplitudes of re- 
sponses to dim light flashes [2|. Fig. 6 shows an example 
of the experimental probability distribution of the am- 
plitudes of the responses to dim flashes of coherent light, 
which clearly indicates the the quantized nature of the re- 
sponses. The dependence was fitted by (4) up to n = 3, 
with Pj (n) considered to be Poissonian with the average 
71 = 1 and the following parameters of P(A\n): Ao = 0.85 
pA,cr A = 0.25 pA, cr D = 0.13 pA. 
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FIG. 6: (color online) Probability distribution of the ampli- 
tudes of the rod response to 80 consecutive flashes of coherent 
light. Saturation amplitude As — 26 pA. Solid black lines are 
given by (3) for n = 0, 1, 2, 3 with the parameters mentioned 
in the text. Dashed blue line is the weighted sum with Pois- 
sonian distribution, given by (4). 

Raw data of the measurements of SNR 

Results of the measurement of SNR for 6 different 
cells, taken from 6 different animals, for the coherent 



and pseudo-thermal light are shown in Fig. 7. Each data 
point is an average response to 80-100 light pulses. 
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FIG. 7: (color online) Dependence of the SNR on the average 
number of impinging photons for coherent (closed symbols, 
solid line) and pseudo-thermal source (open symbols, dashed 
line). Different colors and symbol shapes correspond to dif- 
ferent cells. Each data point is an average response to 80-100 
light pulses. Lines show numerical calculations. 



